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Abstract. Dark-matter halos are the scaffolding around which galaxies and clusters 
are built. They form when the gravitational instability of primordial density fluctu¬ 
ations causes regions which are denser than average to slow their cosmic expansion, 
recollapse, and virialize. Objects as different in size and mass as dwarf spheroidal 
galaxies and galaxy clusters are predicted by the CDM model to have halos with a 
universal, self-similar equilibrium structure whose parameters are determined by the 
halo’s total mass and collapse redshift. These latter two are statistically correlated, 
however, since halos of the same mass form on average at the same epoch, with small- 
mass objects forming first and then merging hierarchically. The structural properties 
of dark-matter dominated halos of different masses, therefore, should reflect this sta¬ 
tistical correlation, an imprint of the statistical properties of the primordial density 
fluctuations which formed them. Current data reveal these correlations, providing a 
fundamental test of the CDM model which probes the shape of the power spectrum of 
primordial density fluctuations and the cosmological background parameters. 


1 The Truncated Isothermal Sphere (TIS) Model 

We have developed an analytical model for the postcollapse equilibrium structure 
of virialized objects which condense out of a cosmological background universe, 
either matter-dominated or flat with a cosmological constant |^,||. The model 
is based upon the assumption that cosmological halos form from the collapse 
and virialization of “top-hat” density perturbations and are spherical, isotropic, 
and isothermal. This leads to a unique, nonsingular TIS, a particular solution 
of the Lane-Emden equation (suitably modified when A ^ Q). The size Vt and 
velocity dispersion cry are unique functions of the mass and redshift of formation 
of the object for a given background universe. Our TIS density profile flattens to 
a constant central value, p^, which is roughly proportional to the critical density 
of the universe at the epoch of collapse, with a small core radius tq ss rt/30 
(where Uy = AttGpoTq and tq = rKing/3, for the “King radius” CKing, defined by 
[Q, p. 228). The density profiles for gas and dark matter are assumed to be the 
same (no bias), with gas temperature T = pmpay/ks- 

These TIS results differ from those of the more familiar approximations in 
which the virialized sphere resulting from a top-hat perturbation is assumed 
to be either the standard uniform sphere (SUS) or else a singular isothermal 
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Table 1 



SUS 

SIS 

TIS|i 

V/Vsvs 

1 

0.833 

1.11;1.07 

T/Tsvs 

1 

3 

2.16;2.19 

Po/pt 

1 

00 

514;530 

{p)lpt 

1 

3 

3.73;3.68 

nlro 

- NA - 

- 00 

29.4;30.04 

Ac/ Ac, SUS 

1 

1.728 0.735;0.774 

K/\W\ 

0.5 

0.75 

0.683;0.690 


The two values refer to flat uni¬ 
verse with A = 0 (left value) and 
flo = 0.3, Ao = 0.7 (right value). 


Fig. 1. (top) Density profile of TIS in a matter-dominated universe. Radius r is in 
units of Tm - the top-hat radius at maximum expansion. Density p is in terms of 
the density psus of the SUS approximation for the virialized, post-collapse top-hat. 
(bottom) Logarithmic slope of density profile. 


sphere (SIS). We summarize their comparison in Fig. 1 and Table I, where rj 
is the final radius of the virialized sphere in units of the top-hat radius at 
maximum expansion (i.e. ?7sus = 0.5), pt = p(d), (p) is the average density 
of the virialized spheres, Ac = (p)/pcrit(^coii), and K/\W\ is the ratio of total 
kinetic (i.e. thermal) to gravitational potential energy of the spheres. 

2 TIS Model vs. Numerical CDM Simulations 

The TIS model reproduces many of the average properties of the halos in nu¬ 
merical CDM simulations quite well, suggesting that it is a useful approximation 
for the halos which result from more realistic initial conditions: 

(1) The TIS mass profile agrees well with the fit to N-body simulations by 
(“NFW”) (i.e. fractional deviation of ~ 20% or less) at all radii outside of 

a few TIS core radii (i.e. outside King radius or so), for NFW concentration 
parameters 4 < cnfw < 7 (Fig. 2). The flat density core of the TIS halo differs 
from the singular cusp of the NFW profile at small radii, but this involves only a 
small fraction of the halo mass, thus not affecting their good agreement outside 
the core. As a result, the TIS central density po can be used to characterize the 
core density of cosmological halos, even if the latter have singular profiles like 
that of NFW, as long as we interpret po, in that case, as an average over the 
innermost region. 

(2) The TIS halo model predicts the internal structure of X-ray clusters 
found by gas-dynamical/N-body simulations of cluster formation in the CDM 
model. Our TIS model predictions, for example, agree astonishingly well with the 
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Fig. 2. Profiles of density (left) and integrated mass (right), for TIS (solid) and NEW 
with cnfw = 4 (short-dashed), 5 (dotted) and 7 (long-dashed) with same (r 2 oo,M 2 oo)- 


mass-temperature and radius-temperature virial relations and integrated mass 
profiles derived empirically from the simulations of cluster formation by P,p^ 
(EMN). Apparently, these simulation results are not sensitive to the discrepancy 
between our prediction of a small, finite density core and the N-body predictions 
of a density cusp for clusters in CDM. Let X be the average overdensity inside 
radius r (in units of the cosmic mean density) X = {p{r))/pb. The radius- 
temperature virial relation is defined as rx = rio{X){T /lO keV)^^^ Mpc, and the 
mass-temperature virial relation by Mx = Mio(X)(T/ 10 keV)^^^M q. 
A comparison between our predictions of rio(A') and the results of EMN is given 
in Fig. 3. EMN obtain Mio(500) = 1.11 ± 0.16 and Mio(200) = 1.45, while our 
TIS solution yields Mio(500) = 1.11 and Mio(200) = 1.55. 

(3) The TIS halo model also successfully reproduces the mass - velocity dis¬ 
persion relation for clusters in CDM N-body simulations and its dependence on 
redshift for different background cosmologies. N-body simulation of the Hubble 
volume [(1000 Mpc)^] by the Virgo Consortium (reported by Evrard at the U. 
of Victoria meeting, August 2000) yields the following empirical relation: 

av = (1080 ± 65) [/i( 0 )M 2 oo/lO^®M 0 ] km/s, (1) 

where M 200 is the mass within a sphere with average density 200 times the 
cosmic mean density, and h{z) = hoy^f2o(l + z)^ -|- Aq is the redshift-dependent 
Hubble constant. Our TIS model predicts: 

av = 1103 [h{z)M2OQ/10^^MQ] km/s. (2) 

(4) The TIS model successfully predicts the average virial ratio, K/\W\, 
of halos in CDM simulations. An equivalent TIS quantity, GM 200 /(r 2 ooav) = 
2.176, is plotted for dwarf galaxy minihalos at z = 9 in Fig. 3(b), from pi , 
showing good agreement between TIS and N-body halo results. A similar plot. 
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Fig. 3. (a) (left) (triangles) Cluster radius-temperature virial relation for CDM simu¬ 
lation results (at 2 ; = 0) as fit by EMN and as predicted by TIS (solid curve). 

(b) GM 2 oo/(uf^’' 20 o) vs. mass for halos from N-body simulations [with Ict error bars]. 
Horizontal line is analytical prediction of TIS model. 


but of Kl\W\ for such halos, was shown by based upon N-body simulations, 
in which the average Rr/|IF| is close to 0.7, as predicted by the TIS model (Table 
I). Those authors were apparently unaware of this TIS prediction since they 
compared their results with the SUS value of Tf/IITl, 0.5, and interpreted the 
discrepancy incorrectly as an indication that their halos were not in equilibrium. 


3 TIS Model vs. Observed Halos 

(1) The TIS profile matches the observed mass profiles of dark-matter-dominated 
dwarf galaxies Q . The observed rotation curves of dwarf galaxies are well fit by 
a density profile with a finite density core given by 

(r/r, + I)(rVr2 + I) 

[||. The TIS model gives a nearly perfect fit to this profile (Fig. 4(a)), with 
best fit parameters Po,b/po,tis = 1-216, Tc/tq^tis = 3.134. This best-fit TIS 
profile correctly predicts Umax, the maximum rotation velocity, and the radius, 
Umax, at which it occurs in the Burkert profile: rmax,B/cmax,Tis = 1-13, and 
Umax,B/Umax, TIS ~ 1.01. 

(2) The TIS halo model can explain the mass profile with a flat density core 
measured by Q for cluster CL 0024-1-1654 at z = 0.39, using the strong grav¬ 
itational lensing of background galaxies by the cluster to infer the cluster mass 
distribution HI- The TIS model not only provides a good fit to the projected 
surface mass density distribution of this cluster within the arcs (Fig. 4b), but also 
predicts the overall mass, and a cluster velocity dispersion in close agreement 
with the value cr„ = 1150 km/s measured by |^. 
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Fig. 4. (a) (left) Rotation Curve Fit. Solid line = Best fit TIS; Dashed line = Burkert 
profile (b) (right) Projected surface density of cluster CL 0024+1654 inferred from 
lensing measurements, together with the best-fit TIS model. 


4 Making Tracks on the Cosmic Virial Plane 

The TIS model yields (poi crv; rt, ro) uniquely as functions of (M, Zcoii)- This de¬ 
fines a “cosmic virial plane” in (po, fo, (Tv)-space and determines halo size, mass, 
and collapse redshift for each point on the plane. In hierarchical clustering models 
like CDM, M is statistically correlated with Zcoii- This determines the distribu¬ 
tion of points on the cosmic virial plane. We can combine the TIS model with 
the Press-Schechter (PS) approximation for Zcoii(M) - typical collapse epoch for 
halo of mass M - to predict correlations of observed halo properties. 

(1) The combined (TIS+PS) model explains the observed Umax “ Umax correla¬ 
tion of dwarf spiral and LSB galaxies, with preference for the currently favoured 
ylCDM model 0 (Fig- 5(a)). 

(2) The TIS+PS model also predicts the correlations of central mass and 
phase-space densities, po and Q = po/cTy, of dark matter halos with their velocity 
dispersions av^ with data for low-redshift dwarf spheroidals to X-ray clusters 
again most consistent with HCDM (Fig. 5 (b,c)). There have been recent 
claims that po =const for all cosmological halos, independent of their mass, as 
expected for certain types of SIDM 10- This claim, however does not seem to 
be supported by most current data (Fig. 5 (c)). 

This work was supported by European Community RTN contract HPRN- 
CT2000-00126 RG29185, grants NASA NAG5-10825 and TARP 3658-0624-1999. 
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mass scale M; typical M{zcou) = M{u = 1)], as labelled with j/-values. All curves are 
(TIS+PS) results, except curve in (a) labelled “Burkert” is a fit to data |^. (a) (left) 
Umax-rmax Correlation. Observed dwarf galaxies (triangles) and LSB galaxies (squares) 
from 0; (b) (middle) Q-ay correlation. Line representing halos of different mass 
which collapse at the same redshift is shown for the case Zcoii = 0. Data points for 
galaxies and clusters are from the following: (1) 49 late-type spirals of type Sc-Im and 
7 dSph galaxies from (open triangles); (^ Local Group dSph Leo I from [0| 

(filled square); (3) 28 nearby clusters, ay from and |^, and po from 0 (crosses). 
(c) (right) Same as (b), except for po vs. ay. 
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